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Abstract: The ^i^SSM is a supersymmetric standard model that solves the /i problem of the 
MSSM using the i?-parity breaking couplings between the right-handed neutrino superfields and the 
Higgses in the superpotential, Xi u^HdHu- The ^ term is generated spontaneously through sneutrino 
vacuum expectation values, ^ = Xi{i)f), once the electroweak symmetry is broken. In addition, the 
couphngs Kijki^fi^ji^i: forbid a global U(l) symmetry avoiding the existence of a Goldstone boson, and 
also contribute to spontaneously generate Majorana masses for neutrinos at the electroweak scale. 
Following this proposal, we have analysed in detail the parameter space of the /iz^SSM. In particular, 
we have studied viable regions avoiding false minima and tachyons, as well as fulfilling the Landau 
pole constraint. We have also computed the associated spectrum, paying special attention to the 
mass of the lightest Higgs. The presence of right and left-handed sneutrino vacuum expectation 
values leads to a peculiar structure for the mass matrices. The most important consequence is that 
neutralinos are mixed with neutrinos, and neutral Higgses with sneutrinos. 
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1. Introduction 



The Minimal Supersymmetric Standard Model (MSSM) is an attractive candidate for 
physics beyond the Standard Model. It not only solves many theoretical puzzles but also 
one expects to find its signatures in the forthcoming large hadron collider (LHC). 

However, in the MSSM lepton and baryon number violating terms in the superpo- 
tential like eab (^XijkLfPjcl + X'-j^LfQ^jdl + fuLfH^^ and X'^jj^d^djUl, respectively, with 
i,j = 1,2,3 generation indices and a,b = 1,2 SU{2) indices, are in principle allowed 
by gauge invariance. As it is well known, to avoid too fast proton decay mediated by the 
exchange of squarks of masses of the order of the electroweak scale, the presence together 
of terms of the type LQd^ and d'^d'^u'^ must be forbidden, unless we impose very stringent 
bounds such as e.g. A'{\2'^ii2 ~ ^ x 10"^"^. Clearly, these values for the couplings are not 
very natural, and for constructing viable supersymmetric (SUSY) models one usually for- 
bids at least one of the operators LQd"^ or u'^d'^d'^. The other type of operators above are 
not so stringently supressed, and therefore still a lot of freedom remains j^]. 

One possibility to avoid the problem of proton decay in the MSSM is to impose im- 
parity conservation (-|-1 for particles and -1 for superpartners). Actually this forbids all 
the four operators above and thus protects the proton. Nevertheless, the choice of -R-parity 
is ad hoc. There are other discrete symmetries, like e.g. baryon triality which only forbids 
the baryon violating operators [^. Obviously, for all these symmetries R-parity is violated. 
Besides, in string constructions the matter superfields can be located in different sectors 
or have different extra U{1) charges, in such a way that some operators violating ii-parity 
can be forbidden but others can be allowed. 

The phenomenology of models where i?-parity is broken differs substantially from that 
of models where i2-parity is conserved. Needless to mention, the LSP is no longer stable, 
and therefore not all SUSY chains must yield missing energy events at colliders. In this 
context the neutralino ^ or the sneutrino Q are no longer candidates for the dark matter 
of the Universe. Nevertheless, other SUSY particles such as the gravitino Q or the axino 
Q can still be used as candidates. Indeed, the well-known axion of the Standard Model 
can also be the cold dark matter. 

There is a large number of works in the literature exploring the possibility of 
i?-parity breaking in SUSY models, and its consequences for the detection of SUSY at 
the LHC [|lO|. For example, a popular model is the so-called Bilinear i?-parity Violation 
(BRpV) model |11|, where the bilinear terms eab lJ'iLfH2 are added to the MSSM. In this 
way it is in principle possible to generate neutrino masses through the mixing with the 
neutralinos without including right-handed neutrinos in the model. One mass is generated 
at tree level, and the other two at one loop. Analyses of mass matrices [^] in the BRpV, 
as well as studies of signals at accelerators have been extensively carried out in the 
literature. Other interesting models are those producing the spontaneous breaking of im- 
parity through the vacuum expectation values (VEVs) of singlet fields ||l^. In the context 



of the Next-to-Minimal Supersymmetric Standard Model (NMSSM) pj, |6|, IT3, |T§ , im- 
parity breaking models have also been studied |19, 20, For a recent review discussing 
the different SUSY models with and without i?-parity proposed in the literature, see ref. 
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There are two strong motivations to consider extensions of the MSSM. On the one 
hand, the fact that neutrino oscillations imply non- vanishing neutrino masses [23|. On 
the other hand, the existence of the /i problem arising from the requirement of a 
SUSY mass term for the Higgs fields in the superpotential, tab l^H^H^, which must be of 
the order of the electroweak scale in order to successfully lead to electroweak symmetry 
breaking (EWSB). In the presence of a Grand Unified Theory (GUT) with a typical scale 
of the order of 10^^ GeV, and/or a gravitational theory at the Planck scale, one should 
explain how to obtain a mass term of the order of the electroweak scale. 



Motivated by the above issues, two of the authors proposed [^, ^g] to supplement 



the superfields i>i contained in the 5f7(2)2,-doublet, Lj, with gauge-singlet neutrino super- 
fields to solve the /x problem of the MSSM. In addition to the usual trilinear Yukawa 
couplings for quarks and charged leptons, and the bilinear /i-term, the right-handed neu- 
trino superfields allow the presence of new terms such as Yukawa couplings for neutrinos 
and possible Majorana mass terms in the superpotential. Besides, trilinear terms break- 
ing i?-parity explicitly such as eah\ i't H^H^ and Kij^OfUji)^ are now also allowed by gauge 
invariance. The // term can be obtained dynamically from the former terms in the superpo- 
tential. When the electroweak symmetry is broken, they generate the // term spontaneously 
through right-handed sneutrino VEVs, // = Aj(P?). In addition, the terms K-ijkt^f'^j^'k forbid 
a global U{1) symmetry in the superpotential, avoiding therefore the existence of a Gold- 
stone boson. Besides, they contribute to generate effective Majorana masses for neutrinos 
at the electroweak scale. Terms of the type v^HdHu and v'^u'^i'^ have also been analysed 
as sources of the observed baryon asymmetry in the Universe |27| and of neutrino masses 
and bilarge mixing ||2^, respectively. 

The superpotencial including the above trilinear couplings with right-handed neutrino 
superfields, in addition to the trilinear Yukawa couplings for quarks and leptons, defines 
the so-called "/x from z/" Supersymmetric Standard Model (;Ui/SSM) |25]. As discussed 
above, the presence of i?-parity breaking couplings in the superpotential is not necessarily 
a problem, and actually the couplings of the //z^SSM are obviously harmless with respect 
to proton decay. Let us also remark that, since they break explicitly lepton number, a 
Goldstone boson (Majoron) does not appear after spontaneous symmetry breaking. As in 
the MSSM or NMSSM, the usual lepton and baryon number violating terms could also in 
principle be added to the superpotential. Actually, even if the terms X-jk^Q'-dt are set 
to zero at the high-energy scale, one-loop corrections will generate them. However, these 
contributions are very small, as we will see in Appendix ^. 

In the ^i/SSM the fi term is absent from the superpotential, as well as Majorana masses 
for neutrinos, and only dimensionless trilinear couplings are present. For this to happen 
we can invoke a symmetry as it is usually done in the NMSSM. Nevertheless, let us 
recall that this is actually what happens in string constructions, where the low-energy 
limit is determined by the massless string modes. Since the massive modes are of the order 
of the string scale, only trilinear couplings are present in the low-energy superpotential. 
String theory seems to be relevant for the unification of interactions, including gravity, and 
therefore this argument in favour of the absence of bare mass terms in the superpotential 
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is robust. 



Let us finally remark that since the superpotential of the ^ui/SSM contains only trilinear 
couplings, it has a Z-^ symmetry, just like the NMSSM. Therefore, one expects to have also 
a cosmological domain wall problem |2£, 3C] in this model. Nevertheless, the usual solution 



1 31] will also work in this case: non-renormalisable operators ||2^ in the superpotential can 
explicitly break the dangerous Z3 symmetry, lifting the degeneracy of the three original 
vacua, and this can be done without introducing hierarchy problems. In addition, these 
operators can be chosen small enough as not to alter the low-energy phenomenology. 

The differences between the /iz^SSM and other models proposed in the literature to 
solve the problem are clear. For example, in the /iz^SSM one solves the problem without 
having to introduce an extra singlet superfield as in the NMSSM, or a special form of the 



Kahler potential |32|, or superpotential couplings to the hidden sector |^, It is also 
worth noticing here that previously studied i?-parity breaking models do not try to address 
the [1 problem. Actually, in the case of the BRpV model the problem is augmented with 
the three new bilinear terms. 

Indeed the breaking of i?-parity generates a peculiar structure for the mass matrices of 
the /iz^SSM. The presence of right and left-handed sneutrino VEVs leads to the mixing of 
the neutral gauginos and Higgsinos (neutralinos) with the right and left-handed neutrinos 
producing a 10x10 matrix. As discussed in ref. |^^, three eigenvalues of this matrix are 
very small, reproducing the experimental results on neutrino masses. Of course, other mass 
matrices are also modified. This is the case for example of the Higgs boson mass matrices, 
where the neutral Higgses are mixed with the sneutrinos. Likewise the charged Higgses 
are mixed with the charged sleptons, and the charged gauginos and Higgsinos (charginos) 
with the charged leptons. 

The purpose of the present work is to extend the analysis of ref. |25], where the 



characteristics of the /iz^SSM were introduced, and only some points concerning its phe- 
nomenology were sketched. Several approximations were considered, and, in particular, 
only one generation of sneutrinos were assumed to acquire VEVs. Here we will work with 
the full three generations, analysing the parameter space of the /ni/SSM in detail, and 
paying special attention to the particle spectrum associated. 

The outline of the paper is as follows: In Section |2| we introduce the model, discussing in 
particular its superpotential and the associated scalar potential. In Section |^ we examine 
the minimisation of the scalar potential. Section ^ is focused on the description of the 
parameter space of the //z^SSM. In Section |5| we thoroughly discuss the strategy followed 
for the analysis of the parameter space and the computation of the spectrum. Section ^ is 
devoted to the presentation of the results. Viable regions of the parameter space avoiding 
false minima and tachyons, as well as fulfilling the Landau pole constraint on the couplings, 
are studied in detail. The associated spectrum is then discussed, paying special attention 
to the mass of the lightest neutral Higgs. Finally, the conclusions are left for Section |^. 
Technical details of the model such as the mass matrices, couplings, one-loop contributions, 
and relevant renormalisation group equations (RGEs), are given in the Appendices. 
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2. The model 

The superpotential of the ^uSSM is given by 

1 



(2.1) 



where we take HJ = (Hj, H^), = {H+, H^), Qj = {in, di), Lf = e,), and Y, A, k, 
are dimensionless matrices, a vector, and a totally symmetric tensor, respectively. In the 
following the summation convention is implied on repeated indices. 

In order to discuss the phenomenology of the /ii^SSM, and working in the framework 
of gravity mediated SUSY breaking, we write the soft terms appearing in the Lagrangian, 
-^soft) as 



-£soft = ry?Q^^ Qf Q'^ + uf u] + df d'] + m| Lf + m^c. ef 

2 Tja* Tja I _2 zja* rja , _2 



J 



(AM,, Hi Q1 u] + (A,y,),, H'i Q\ d] + {A,Y,\, L\ ~e] 



+ (^,y,),,F^L^p| + c.c. 



+ 



ki'tv^M +C.C. 



M3 A3 A3 + M2 A2 A2 + Ml Ai Ai + c.c. . 



(2.2) 



In addition to terms from £soft, the tree- level scalar potential receives the usual D and 
F term contributions. Thus, the tree-level neutral scalar potential is given by 



(2.3) 



where 



+ [av^.^lhv] - ax^vlHlHl + ^a«^^,Pfz>|P^ + c.c. 



rO u-O* 



2 ~c,~c* 



with Ci.^^. = (A,,i;.)ij, = {Ax\)i, a^^^^ = {A^K)ijk, 



(2.4) 



(2.5) 



with G'^ = gf + §2, and 

- {'^ijk>^*jHa*H^*i>fi>k - YuijK'ljkH^i>ii'f*i'k* + Yu^^ X*j Hj* H^* H^i>i 

+ y4A,FK^>*z>r + c.c.) 
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Once the electroweak symmetry is spontaneously broken, the neutral scalars develop 
in general the following VEVs: 



{H^,) = v,, {Hl)=vu, {h) = v^. ^ = ^. (2.7) 



In the following we will assume for simplicity that all parameters in the potential are real. 
Although in 'multi-Higgs' models with real parameters the VEVs of the neutral scalar fields 
can be in general complex [^], the analysis of this possibility is beyond the scope of this 
work, and we leave it for a forthcoming publication, where spontaneous CP violation will 



be studied in detail [36|. Nevertheless, it is worth noticing here that this assumption of real 
VEVs is consistent once one guarantees that the minimum with null phases is the global 
one. It is straightforward to see that this is guaranteed in general for the VEVs v^, fd, 
, imposing the conditions Aj > 0, Km > 0, Ax- > 0, A^^.^ < 0, and ^^^^.^ = K^-fc = if 
i = j = k is not satisfied. Concerning the VEVs fj, it is sufficient to impose Y^.- > 0, and 
Y,^.^ = Au^. = for z 7^ J, with the extra condition 

+ XjU^jU^Vd - A^^^Vui^j - KijkiyjUkVu > . (2.8) 



The above conditions on the signs of the parameters, together with ( |2.8| ), will be used for 
the analysis of the parameter space and spectrum of the /ii^SSM in Section 0. 



3. Minimisation of the potential 



As mentioned above, the EWSB generates the VEVs written in eq. (p. 7]). Thus one can 
define as usual 



H^ = hu+ iPu + Vu , H^ = hd+ iPd + Vd, 
= {^t)'' + ii^tY + ^1 , n = i^i)'' + ^i^i^^Y + n ■ (3.1) 

Then, the tree-level scalar potential contains the following linear terms: 

<ear = t>d + tK + (?f )« + 1%^ (z7,)« , (3.2) 

where the different ^ are the tadpoles at tree-level. They are equal to zero at the minimum 
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of the tree-level potential, and are given by 

+ XiXiVdvi - XjKijkVuvfh'l - Yi^^^Xki^iU^i^j - Yy^^Xjvlvi , (3.3) 

+ XiXjVul^iiyj + XjXjvjvu - XjKijkVdVil^k + ^v.jK'ljkVil'tl^l 

- 2XjYy^^VdVuVi + Yy^Yu^^Vuvlv^j + Yy^Yi,^^VuyiVk , (3.4) 
tic = m~c^i^j + ay.MjVu - ax^VuVd + a^^^.^Ujuf. + XiXjvluj + XiXjvjuj 

- 2XjKijkVdVuVl + 2KiimKljkiy^UjV^ - Y^^^Xkl^jU^Vd - Y^^.XiVdl^kl^j 

+ 2Y^^^KikiVui^ji'f + Y^^Yy^^VjUivl + Y^^Y^^^vlv^j , (3.5) 

t% = ^G^i^ji'j + vl)vi + m^^^Uj + Qy^^Vuv'-j - Y^^^XkVdV^jvl 

- Yu,,Xjvlvd + Y^^^KijkVui^ji^k + Yu^Y^^i^i^iiy^i^^ + Y^^^Y^^^vlvj . (3.6) 

As it is well known, in order to find reliable results for the EWSB, it is necessary to 
include the one-loop radiative corrections. The effective scalar potential at one-loop level 
is then 

V = V° + V^, (3.7) 
where includes bubble diagrams at one-loop with all kinds of (s)particles running in the 



loop |37]. Minimizing the full potential is equivalent to the requirement that the one-loop 
corrected tadpoles, t = + 1^, where represents the one-loop part, vanish. 

Let us finally remark that, since minima with some or all of the VEVs in eq. ( |2.7] ) 
vanishing are in principle possible, one has to check that the minumum breaking the elec- 
troweak symmetry, and generating the /i term spontaneously, is the global one. This will 



be studied in detail when analyzing the parameter space of the model in Section 6.1 
4. n^SSNl parameter space 

At low energy the free parameters in the neutral scalar sector are: Aj, Kijk, fnu^, f^Hui 
f^L-i "^^'fj' -^^^■> -^i^rjk-: ^'^d A^ij- Strong upper bounds upon the intergenerational scalar 
mixing exist |38], so in the following we assume that such mixings are negligible, and 



therefore the sfermion soft mass matrices are diagonal in the flavour space. This occurs 
for example in several string compactifications as a consequence of having diagonal Kahler 
metrics, or when the dilaton is the source of SUSY breaking |^9|. Thus using the eight 
minimization conditions for the neutral scalar potential in the previous section, one can 
eliminate the soft masses mn^-, rriH^, mj^ , and m^rc in favour of the VEVs Vd, Vu, Vi, and 
f?. On the other hand, using the Standard Model Higgs VEV, v 174 GeV, tan/3, and 
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Ui, one can determine the SUSY Higgs VEVs, Vd and Vu, through = + + ^i- We 
thus consider as independent parameters the following set of variables: 

f^ijki tan/3, (4.1) 

It is worth remarking here that the VEVs of the left-handed sneutrinos, fj, are in 
general small. Notice that in eq. ( |3.6[ ) ^ as — > to fulfil t^. = 0, and since the 
couplings Yi, determine the Dirac masses for the neutrinos, Y^Vu ^ rriD < 10~^ GeV, the 
z^'s have to be very small. Using this rough argument one can also get an estimate of the 



values, < m£) |25|. Then, since Vi « Vd,Vu we can define the above value of tan/3 as 
usual, tan (3 = —. 

Assuming for simplicity that there is no intergenerational mixing in the parameters of 
the model, and that they have the same values for the three families (with the exception 
of for which we need at least two generations with different VEVs in order to guarantee 
the correct hierarchy of neutrino masses), the low-energy free parameters in our analysis 
will be 

A, K, tan/?, vi, 1^3, u", Ax, A^, Ay , (4.2) 



where we have chosen v\ = V2 v^, and we have defined A = Aj, k = Km, v'^ = uf, 
Ax = Axi , A^ = Ak^.. , Ay = Ay-- . Nevertheless, let us remark that the formulas given in 
the Appendices are for the general case, without assuming universality of the parameters 
or vanishing intergenerational mixing. 

The soft SUSY-breaking terms, namely gaugino masses, Mi^2,3! scalar masses, mg g 
and trilinear parameters, Au^d,ei ^-re also taken as free parameters and specified at low scale. 
Data on neutrino masses, and the usual Standard Model parameters such as fermion and 
gauge boson masses, the fine structure constant a{Mz), the Fermi constant from muon 



decay G^, and the strong coupling constant as{Mz), will be used in the computation |40|. 
Concerning the top mass, we will take mt = 172.6 GeV [^]. 

5. Strategy for the analysis 

We now show the algorithm used in the analysis of the model. In particular in the analysis 
of the parameter space, and in the computation of the spectrum. Below Mz, a{Mz) and 
as{Mz) are first evolved to 1 GeV using 3 loop QCD and 1 loop QED Standard Model 



/^-functions |42]. Then the two gauge couplings and all Standard Model fermion masses 
except the top quark mass are run to Mz- The /^-functions of fermion masses are taken to 
be zero at renormalisation scales below their running masses. The parameters at Mz are 
used as the low energy boundary condition in the rest of the evolution. 

We work in the dimensional reduction (DR) scheme |43| in which the counterterms 



cancel only the divergent pieces of the self-energies required to obtain the pole masses. 
Thus, they become finite depending on an arbitrary scale Q and the tree level masses are 
promoted to running masses in order to cancel the explicit scale dependence of the self- 
energies. It implies that all the parameters entering in the tree-level masses (couplings and 
soft masses) are DR running quantities. 
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The algorithm proceeds via the iterative method, and therefore an approximate initial 
guess of the /iz^SSM parameters is required. As explained above, from tan/? and Mz one 
can determine the Higgs VEVs, Vd and Vu, and from these the third family DR Yukawa 
couplings can be approximated as 

= ^ , nAQ) = , (5.1) 

Vu Vd 



where Q = mt{mt) is the renormalisation scale. The MS values of fermion masses are used 
for this initial estimate. The fermion masses and Og at the top mass scale are obtained 
by evolving the previously obtained fermion masses and gauge couplings from Mz to mt 
(with the same accuracy). The electroweak gauge couplings are estimated by ai{Mz) = 
5a{Mz) /S cos'^ 6w, ct2(-^z) = a{Mz)/ sin^ 9w- Here, sin^^y is taken to be the on-shell 
value. These two gauge couplings are then evolved to with one-loop Standard Model 
/3-functions. 

The gauge and Yukawa couplings and the VEVs are then evolved to the scale (in the 
first iteration we guess Ms) 



Ms ^ ^mi^{Ms)mi^{Ms) , (5.2) 
where the scale dependence of the electroweak breaking conditions is smallest |E3]. For 



it we employ the one-loop DR /3-functions given in Appendix E. The supplied boundary 
conditions on the soft terms are then applied. 

At this point we determine the neutrino Yukawa couplings through the 10 x 10 neutral 



fermion mass matrix which can be written as |25| 




Mn = \ r n ' (5-3) 



where M is a 7 x 7 matrix composed by the MSSM neutralino mass matrix and its mixing 
with the f ?, while m is a 7x 3 matrix containing the mixing of the vi with MSSM neutralinos 



and the f?. The full matrix is written in Appendix A. 3. 

The above matrix is of the see-saw type giving rise to the neutrino masses which in 
order to account for the atmospheric neutrino anomaly have to be very small. This is the 
case since the entries of the matrix M are much larger than the ones in the matrix m. 
Notice in this respect that the entries of M are of the order of the electroweak scale while 
the ones in m are of the order of the Dirac masses for the neutrinos. Therefore in a first 
approximation the effective neutrino mixing mass matrix can be written as 

meff = — m"^ • M~^ ■ m . (5-4) 

Because mefr is symmetric and mggmefr is Hermitian, one can diagonalise them by a unitary 
transformation 

t^MNS "T-eff f^MNS = diag {m^^ , 1X1^., , m^s) , (5.5) 
^MNS "iefr"^eff f^MNS = diag (mu^,m^^,m^.^) . (5.6) 
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The masses are connected with experimental measurements through 



= \l + ^ml^i , = \J ml, + Am^^^ . (5.7) 

To determine the neutrino Yukawa couphngs we choose the basis where 1^ is diagonal. 
Then we employ a numerical procedure which consists in solving three non-linear coupled 
equations in Yj^. - determined by the diagonalisation of me/ / . Another way would consist in 
fixing the neutrino Yukawa couplings as inputs giving the left-handed sneutrino VEVs as 
outputs. However the method employed is appropriated from the numerical stability point 
of view. 

The determination of the charged lepton Yukawa couplings should follow a similar 



procedure through the charged fermion mass matrix written in Appendix A. 2. In this 
matrix the charginos are mixed with the charged leptons. However, because it turns out 
that Y^.. < 10~^ in order to achieve the smallness of the neutrino masses (and also < 10~^ 
GeV as discussed in Section the 2x2 chargino submatrix is basically decoupled from 
the 3x3 charged lepton submatrix. Thus the charged lepton Yukawas can be determined 



directly from the charged lepton masses |]2^, as it is stated above in eq. (5.1). 

At the Ms scale the tree-level tadpoles, eqs. (3.4-3.7), are set to be zero to guarantee 
the EWSB. As discussed above, at this scale the scale dependence of the EWSB parameters 
is smallest. In this way the soft masses {Ms), mjj (Ms), m- {Ms), and m~c{Ms), are 
derived. 

The next step consists of performing a check in order to ensure that the minimum which 
breaks the electroweak symmetry spontaneously is the global one. For it, we compute the 



global minimun using a 'genetic' algorithm for global optimisation [45| which has a high 
performance. Then we compare it with the physical one. 

In the final step the DR (tree- level) superparticle mass spectrum consisting of squarks, 
CP-even (odd) neutral scalars, charged scalars, neutral fermions and charged fermions (see 
Appendix A) is determined at the Ms scale. Notice that once the tree-level mass spectrum 
is known, radiative corrections to the neutral scalar potential and the tadpoles as well as 
for computing pole masses are calculable. 

In order to check the absence of a Landau singularity (by requiring any Yukawa cou- 
pling to be less than \/4 tt) the Yukawa couplings are evolved to the GUT scale. Finally, 
Yukawa couplings, gauge couplings, and VEVs, are evolved back down to Mz, and SUSY 
one- loop thresholds containing squark/gluino in the loop are added to the third family of 
quark Yukawa couplings and to the strong coupling constant [p6|| . The whole process is 
iterated, as it is sketched in Fig. ||, with the inclusion of one-loop corrections to the neutral 
scalar potential. It is equivalent to add the one-loop tadpoles to eqs. (3.4-3.7). Then 
the global minimun is computed following the procedure described above. For this work 
we have computed the leading one- loop contributions to the tadpoles, which come from 
(s)quarks in the loops, in the DR scheme. The results are given in Appendix C. For the 
neutral scalar potential we employ the results in ref. |37|. 



Once the DR sparticle masses all converge to better than the desired fractional ac- 
curacy, the computation of the physical masses requires the addition of loop corrections. 
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Set VEVs, Yukawas, and add SUSY rad. corr. to gs{Mz), ht,b{Mz) 



Run to Ms 



Apply soft SUSY-breaking boundary conditions 



Neutrino Yukawas, EWSB 



Run Yukawas to Mqut- Check for Landau Poles 



Run to Mz 



Run to Mg. Calculate sparticle pole masses 



Figure 1: Iterative algorithm used to calculate the SUSY spectrum. Each step (represented by a 
box) is detailed in the text. The initial step is the uppermost one. Ms is the scale at which the 
EWSB conditions are imposed, as discussed in the text. 



It is well known that the role of the radiative corrections to the lightest CP-even Higgs 
boson mass is extremely important (see ref. for studies of this effect in the NMSSM). 
The leading ones come from an incomplete cancellation of the quark and squark loops. 
Following this we have added those corrections as described below. The rest of the masses 
are tree-level DR running masses. 
The gluino mass is then given by 

mf^ = MsiMs) . (5.8) 

The rest of SUSY particles mix in the interaction basis, and a rotation to their mass 
states basis is required. The scalar sector includes the squarks which are MSSM-like and 
therefore their masses and mixing angles are the result of performing a Jacobi 2x2 rotation 



of the matrices in Appendix A. 1.4. Charged and neutral fermion masses are the result of 
diagonalising their mass matrices which are given in Appendices |A.2| and [A.3| respectively. 

It is worth mentioning that a final check is required to see if the procedure used above 
to compute neutrino Yukawa couplings is consistent with the final lightest eigenvalues of 
the neutralino mass matrix eq. ( |5.3| ). 

The CP-even scalar masses are obtained from the real parts of the poles of the propa- 
gator matrix 

Det [pjl - Mso (p?)] = , mj^ IZeipf) , (5.9) 
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where 

Msoip'') = M^iQ) + Usoip^ Q) , (5.10) 

with n'^° being the matrix of the renormahsed self-energies in the DR scheme of the CP- 
even scalars. The ones involving quarks and squarks in the loop are shown in Appendix D. 
We diagonalise the matrix .Msoipf) at an external momentum scale equal to its pole mass 
pf = mf through an iterative procedure. 

Finally, the quark Yukawa couplings, gauge couplings, and VEVs are evolved back 
down to Mz- 



6. Results and discussion 

Using the results of the previous Sections and Appendices, we will study in detail the 
parameter space and spectrum of the nvSSM. 

6.1 Analysis of the parameter space 

In this subsection the parameter space of the fiuSSM will be studied. We will see that 
avoiding the existence of false minima and tachyons, as well as imposing perturbativity 
(Landau pole condition) on the couplings of the model, important constraints on the pa- 
rameter space will be found. 

The free parameters of our model have already been presented in eq. ( |4.2| ). As afore- 
mentioned, we take them to be free at the electroweak scale. As discussed in Section ^, we 
will determine the neutrino Yukawa couplings through the experimental data on neutrino 
masses. We will use the direct hierarchical difference of masses, taking the typical values 
TTij,^ = 10"^2 QgY^ = g i X 10-12 GeV and 7711,3 = 4.7 x 10"^^ GeV. Finally, as discussed 
in Section |^, it is sufficient to work with only two different left-handed sneutrino VEVs. 
In particular, we choose vi = 1^2 = 1-4 x 10~^ GeV and = 1.4 x 10~^ GeV, which are 



typical values in order to satisfy the minimum equations (3^) and data on neutrino masses 
through the see-saw mechanism ( [5. 3D . Possible variations of these values will not modify 
qualitatively our results below. 

Throughout this section we will consider several choices for the values of 

A, K, tan/?, i^'', Ax, A,^ ,Au , (6.1) 

using the sign conditions explained in Section ^. Besides, we work with a negative value of 
j4jy in order to fulfill condition ( |2.8D more easily. 

Concerning the rest of the soft parameters we will take for simplicity in the computation 
rriQ -c gc = 1 TeV, Au^d,e = 1 TeV, and for the gaugino masses only M2 = 1 TeV will be 
used as input, whereas the others will be determined by the approximate GUT relations 
Ml = ^M2, M3 = ^Ms, implying Mi O.SMa, M3 2.7M2. 

Let us first discuss when the minimum we find following Sections § and |^ is the global 
one. In particular, one has to be sure that it is deeper than the local minima with some 
or all of the VEVs in eq. (|2.7| ) vanishing. Concerning the latter one can check that the 
most relevant minima are the solutions with only Vu or i/'^ different from zero (in some 
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special situations also the case with all VEVs vanishing can be relevant). For example, for 
a given value of v'^ the term proportional to in (2.4) turns out to be important: the 
more negative the value of A^, the deeper the minimum becomes. This might in principle 
give rise to a value of the potential ( |2.3D in the direction with only i/^ 7^ 0, more negative 
than the one produced in the realistic direction with all VEVs non vanishing. In that case 
the associated points in the parameter space would be excluded by the existence of false 
minima. Notice that ni^j^ is independent on the value of as can be deduced from eq. 
(3.4) with = 0. Thus although ni^u^ will contribute to the realistic direction, it plays no 
role in the above argument. 

On the other hand, we can also deduce from eq. (3.4) that for reasonable values of the 
parameters the larger u^, the smaller rn^fj^ becomes in order to cancel t^. As a consequence, 
the realistic direction becomes deeper, and the associated points in the parameter space 
are allowed. 

Both effects can be seen in Fig. ||a, where the (^k, i^^) parameter space (recall our 
assumption = v^) is plotted for an example with A = 0.1, k = 0.4, tan/3 = 5, and 
Ax = —Ay = 1 TeV. For a given value of we see that for A,^ sufficiently large and 
negative one obtains a false minimum (gray area). For larger values of v'^ one needs values 
more negative of A^^ to obtain the false minimum. Let us remark that although m?c depend 
on as can be obtained from eq. (3.5), we can see in Fig. |2|b that this variation is not 
crucial for the discussion above. Notice that the values of m?c for points of the parameter 
space close to the false minimum area do not vary in a relevant way. 

In Fig. I we can also see that part of the parameter space is excluded due to the 
occurrence of tachyons in the CP-even neutral scalar sector. Thus the realistic direction 
with all VEVs non-vanishing is not even a local mimimum. This happens in general 
when the off-diagonal values \M? .-cnrI or |M? .-cnbI of the CP-even neutral scalar matrix 
(see App endix [A . 1 . l| ) become significantly larger than |M^~c')fl(-~c-)fl| in some regions of the 
parameter space, thus leading to the appearance of a negative eigenvalue. The violet area 
in Fig. |2| corresponds to this situation. In particular, notice that the relevant terms in the 
off-diagonal pieces are linear in z/'^, whereas in M^-c-jfl^-c-jjj they are quadratical. Thus, for 
a given value of yl^; the smaller the value of v'^, the smaller the latter terms become giving 
rise to the possibility of tachyons. Notice also that there is a term proportional to in 
M^~c implying that, for a given value of v"^, the more negative the value of A,^, the 

smaller M^-^^h^-c-^r become. This is also reflected in Fig. ^. 

Let us now discuss the possibility of minima deeper than the realistic one in the 
direction with only Vu 7^ 0. When the values of v'^ are large, we can see from eq. (3.5) 
that m?c must be negative in order to cancel the cubic term in u'^. However, when the 
values of v"^ are small, m|c must be positive in order to cancel the quadratic term in v'^ 
proportional to a^, which is now the relevant one. This may give rise for small i''^ to a 
value of the potential ( |2.3| ) in the direction with only Vu 7^ 0, more negative than the one 
produced in the realistic direction with all VEVs non vanishing. This situation is shown in 
Fig. ^ where the {Ax, y^) parameter space is plotted for an example with A = 0.1, k = 0.4, 
tan/3 = 5, A^^ = Ay = —1 TeV. We can see in Fig. ^3 that the smaller v'^, the larger r?T,?c 
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(a) (b) 

Figure 2: (A^, i^'^) parameter space for tan/3 = 5, A = 0.1, k = 0.4, and Ax = —A^, = 1 TeV. In both 
cases the gray and violet areas represent points which are excluded by the existence of false minima and 
tachyons, respectively. In (a) the colours indicate different values of the soft mass m-l^ . In (b) the 
colours indicate different values of the soft masses mL. 



become, making it easy the appearance of a false minimum. Let us remark that the points 
in the gray area above ~ 1 TeV are actually forbidden by minima deeper than the 
realistic one with all VEVs vanishing. Notice to this respect in the figure that those points 
correspond to positive values of m'jj^ and m?c. This is also true for Fig. ^ discussed below, 
but for points above Ax ^ 2 TeV. 

It is worth noticing here that m?c is essentially independent on the value of Ax, as can 
be easily deduced from eq. (3.5). On the other hand, we can see from eq. (3.4) that rrijj^ 
does depend on Ax through the term proportional to ax- In particular, if we decrease Ax, 
rn^j^ also decreases, as shown in Fig. Although this might in principle contribute to 
produce a minimum deeper than the realistic one in the direction with only t>„ 7^ 0, we 
see in the figure that for the parameter space studied the effect is negligible. Nevertheless, 
increasing the value of A, ax also increases, and this effect can be more important. This is 
shown in Fig. where A = 0.2 is considered. We can see that the parameter space is now 
more constrained. We also show in Fig. ^ the values of m?c in the allowed region. 

Actually, there is a new tachyonic region for large values of v'^. This happens because 
the off-diagonal value in Appendix ( |A.l.l| ) has a quadratic dependence on v'^, thus 

leading to the appearance of a negative eigenvalue. Notice in this respect that a similar 
dependence in the diagonal pieces \^h^hj\ ^"^^ \-^huhu\ canceled once we substitute the 
value of the soft masses using eqs. (3.3) and (3.4). 
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Figure 3: {A\, v"^) parameter space for tan/3 = 5, A = 0.1, k = 0.4, and = Ay = —1 TeV. In both 
cases the gray and violet areas represent points which are excluded by the existence of false minima and 
tachyons, respectively. In (a) the colours indicate different values of the soft mass m^^ . In (b) the 
colours indicate different values of the soft masses mic. 
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Figure 4: The same as in Fig. 3 but for A = 0.2. 
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For each point in the parameter space, one also requires perturbativity, i.e. the absence 
of Landau singularities for the couplings. Let us discuss now in detail the case of A, since 
this is the relevant coupling when discussing the upper bound on the lightest Higgs mass 
in the next Subsection. 

Once perturbativity is imposed, the value of A is actually bounded. To obtain a rough 



estimation we can use eq. (E.ll) in the Appendix neglecting Y^. ., and taking Km = k and 



K-ijk = if i = j = A: is not satisfied. Then we can write that equation as 

where we have defined = AjAj, i = l,...,n, with n the number of singlets, and C is 
a quantity independent on Aj. It is worth noticing here that the RGE for the relevant 
parameter is clearly independent on n. Thus we could in principle expect a bound for 
A'^ similar to the one of the NMSSM for A. Recall that in the NMSSM there is only one 
singlet, and A^ < (0.7)^. To complete the discussion we can solve a simplified version of 
eq. ( |6.2D neglecting the piece proportional to C, with the result 

A^W) = , il^f^n,. > (6.3) 
1^ 



1 + 



where Q is the renormalization scale, and Qq the scale of the high-energy theory. At 
the high-energy scale the Landau pole condition for each coupling can be imposed as 
^j{Qo) < 47r, implying A^(Qo) < 47rn, and therefore one obtains the following upper 
bound: 

4n7r 

1 + 



< ^^^TTToiTT ■ (6-4) 



For Qq sufficiently large the second term in the denominator is much larger than one, and 
the equation can be approximated as 

For example, if the high-energy theory is a typical GUT with Qq ~ 10^^ GeV, then from 
eq. ( |6.5| ) with Q ~ 100 GeV one obtains the low-energy bound A^ < (0.78)2. taking into 
account that C in eq. ( |6^ ) gets a negative (positive) contribution from the top (gauge) cou- 
pling, one should expect a final bound slightly stronger. The numerical analysis indicates 
that this is the case, with A^ < (0.7)^ as expected. Thus in our case where i = 1,2, 3, we 
obtain the bound for each coupling A = Aj < 0.7/\/3 ^ 0.4. 

Although in the numerical analysis below we will impose the Landau pole constraint 
assuming that the perturbative description of the model is valid up to the GUT scale, it 
is worth noticing here that intermediate scales like 10^^ GeV seem also to be interesting 
to explain several experimental observations. In addition, it has been found that the 



string scale may be anywhere between the weak and the Planck scale |49]. Also NMSSM 



like models restricted to be perturbative up to about 10-100 TeV have been studied pO| 
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Considering these possible uncertainties in the unification scale, and using e.g. Qo ~ 10 
GeV, from eq. ( |6.5D we would obtain A'^ < (0.95)^. Taking into account as above the 
other contributions to the RGE, one can find the final bound < (0.88)^, and therefore 
Xi < 0.88/\/3 ~ 0.5. It is worth noticing then that, for intermediate scales the allowed 
parameter space is larger than in the case of a typical GUT. Obviously, smaller scales 
would imply even larger allowed regions. For example, with Qo ~ 10 TeV, one obtains a 
final bound A'^ < (1.91)^, implying Aj < 1.1. Another modification will be related to the 
lightest Higgs mass. As will be discussed in the next Subsection, its upper bound is also 
larger for smaller unification scales. 

In Figs. |5[]^ we study the (A, k) parameter space. As expected from the above dis- 
cussion, A < 0.4. Concerning the value of k, we also see that perturbativity up to the 
GUT scale imposes the bound k < 0.6, similarly to the NMSSM. In Fig. |5| we show an 
example with tan P = 5, Ax = —A^ = —Ai, = 1 TeV, and = 2 TeV. For A > 0.05 a false 
minimum region appears. As we can deduce from Fig. |5|a, the reason is that m'jj^ becomes 
large and negative, producing as a consequence a minimum deeper than the realistic one 
in the direction with only Vu 0. 

It is clear from Fig. ^ that the presence of tachyons increases for large values of A (see 
e.g. the orange area). The reason is that the off-diagonal value |M^^^^| in Appendix ( |A.l.l| ) 
has a dependence on ax, thus leading to the appearance of a negative eigenvalue. We can 
also see to the left of the figure, for very small values of A, a narrow band with tachyons. 
The relevant off-diagonal piece is now {u'=)r\' Notice that there are terms with opposite 
signs producing a cancellation of the mixing for particular values of A. However for very 
small values the cancellation disappears and a large mixing producing negative eigenvalues 
arises. 

In Fig. ^ we show the modifications produced by a decrease in the value of v'^. In 
particular, we consider the same values of the parameters as in Fig. ^ but with u'^ = 1 TeV 
instead of 2 TeV. The allowed region is now reduced. Notice that m?c becomes positive for 
larger values of k, producing the presence of minima deeper than the realistic one in the 
direction with only Vu 7^ 0. Let us also remark here that the points in the gray area about 
A ~ 0.05 and k ~ 0.35 are actually forbidden by minima deeper than the realistic one with 
all VEVs vanishing. 

Decreasing further the allowed region decreases, and in particular for z^^ « 500 GeV, 
and the same values of the parameters as above, we find that the whole region disappears. 
Nevertheless this situation can be improved if we modify the values of A^ and Ax- In 
particular, decreasing Ax, and increasing (decreasing in modulus) Af^, the terms in the 
potential proportional to them contribute to generate a realistic minimum. This can be 
seen in Fig. |a, where we take = 500 GeV, = 200 GeV, and yl«; = -50 GeV. The 
allowed region is even larger than in Fig. ^ where u'^ = 1000 GeV. 

Let us finally discuss the variation in tan /3. Larger values of tan (3 lead to an increase 
of the mixing in the CP-even neutral scalar matrix, and as a consequence the tachyonic 
region is larger. We show this effect in Fig. ^b for tan (3 = 20. Although the allowed 
region is smaller than in Fig. |6|, the effect is not very important. This is also true for larger 
values of tan/3. The reason being that the large value of i^'^ = 1 TeV produces a heavy 
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(a) (b) 

Figure 5: (A, k) parameter space for tan/3 = 5, A\ = 1 TeV, = A,^ ~ —1 TeV, and = 2 TeV. In 
both cases the gray and violet areas represent points which are excluded by the existence of false minima 
and tachyons, respectively. The yellow area represents points which are excluded due to the occurrence 
of a Landau pole. The orange area is excluded by both, Landau pole and tachyons. In (a) the colours 
indicate different values of the soft mass . In (b) the colours indicate different values of the soft 
masses mL . 




(a) (b) 
Figure 6: The same as in Fig. 5 but for the case v"^ = 1 TeV. 
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Figure 7: The same as in Fig. 5 but for the cases (a) tan/3 = 5, Ax = 200 GeV, = —50 GeV, 
Ai, = —1000 GeV, and ly'^ = 500 GeV. The colours indicate different values of the soft mass rnj^ . 
(b) tan/3 = 20, Ax = 1000 GeV, A,, = A^ = -1000 GeV, and ly" = 1000 GeV. The colours indicate 
different values of the soft masses toL. 



right-handed sneutrino, and therefore a large entry |M^~c^^^~c^^|. Since the other relevant 

entries, iM^hl and \M? /~c-,ii|, are generically much smaller, it turns out to be difficult to 
generate a negative eigenvalue. As for tan /3 = 5, decreasing further i^^ for the same value of 
the parameters, the allowed region decreases. Both effects, the generation of false minima 
and tachyons, are contributing significantly to forbid points of the parameter space. In 
particular, the latter effect also contributed to forbid the whole region for tan /3 = 5 and 
v'^ ~ 500 GeV. This is obvious, since the potential is bounded from below, and, as a 
consequence, the existence of tachyons implies the existence of a deeper minimum. The 
whole region is also fobidden for tan P larger than 5 when i/^ ~ 500 GeV. 

6.2 Analysis of the spectrum 

Let us now discuss general characteristics of the particle spectrum of the /ii^SSM. The 
breaking of i?-parity generates a peculiar structure for the mass matrices. The presence 
of right and left-handed sneutrino VEVs leads to mixing of the neutral Higgses with the 
sneutrinos producing the 8x8 neutral scalar mass matrices for the CP-even and CP-odd 
states written in eqs. ( |A.l ) and ( |A.14 ), respectively. Note that after rotating away the 



CP-odd would be Goldstone boson, we are left with seven states. It is also worth noticing 
here that the 5x5 Higgs-right handed sneutrino submatrix is basically decoupled from 
the 3x3 left handed sneutrino submatrix, since the mixing occurs only through terms 
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proportional to or Y^y^^, which are therefore neghgible. 

Given the interest of the hghtest Higgs boson mass in the analysis of SUSY models, it 
is worth discussing here its upper bound in the iivSSM.. Let us recall that for an extension 
of the MSSM with singlets Si, i = generating the fi term through the couplings 

fiafe-^i Si H^H^, one can obtain a tree-level upper bound on the lightest neutral Higgs mass 
1 51, 52 1 using the 2x2 sub matrix defined by ruHu m//^ (see Appendix [A.l.l| ), 



ml < M| (^cos2 2P + sin2 2(3^ « M| (cos^ 2/3 + 3.62 sin^ 2/3) , (6.6) 

where = AjAj was defined in the previous Subsection. Neglecting the small neutrino 
Yukawa couplings Fj^-^ and with the substitutions Si — > u^, i = 1,2, 3, the superpotential of 
the /ii^SSM (|2.lD is equivalent to the above extension, and therefore we can use the same 
bound (|6.6|) in our computation. 

Clearly, one can optimise this bound choosing tan (3 as small as possible, as well as 
A as large as possible. Concerning the latter, let us recall our discussion in the previous 
Subsection: the value of A is actually bounded once perturbativity is imposed, and, in 
particular, we found A^ < (0.7)^ for a typical GUT. Now, using this bound one can write 
(|6^) as 



I < M| (cos^ 2/3 + 1.77 sin^ 2/3) , (6.7) 

which indicates that for small values of tan (3 (i.e. large values of sin 2/3) one might obtain 
in principle large tree-level values for the lightest Higgs mass, unlike the MSSM where the 
second term in ( |6.7D is absent. For example, for tan /3 = 2(4) one obtains m/j < 1.22(1.08) x 
Mz ~ 111(98) GeV. 

Of course, in order to get masses close to the upper bound, choosing a certain range of 
values for other parameters of the model in (^]^) is also necessary. In particular, we must 
avoid as much as possible the mixing of the light eigenstate h of the 2x2 Higgs submatrix 
in Appendix ( |A.l.lD with the right-handed sneutrinos (see eqs. (IXsl) and ([Zel) ). bmce 
this submatrix is essentially diagonalized by the angle ^ — /?, it is easy to check that one 
has to impose 

XlGXiy" - {Ax + 2ki/^) sin 2/?] ^ . (6.8) 
On the other hand, it is well known that the one-loop correction to the lightest Higgs 



mass can be very important. One can check that, similarly to the NMSSM |48], the upper 
bound for the lightest doublet-like Higgs mass of the /iz^SSM is of the order of 140 GeV for 
tan /3 ~ 2. 

As discussed in the previous Subsection, for high-energy theories with smaller funda- 
mental scales than the GUT one, the upper bound for the coupling turns out to be larger. 
In particular, for an intermediate scale of the order of 10^^ GeV we found A^ < (0.88)^. 
Thus, from (|6.6|), one is also able to get a larger tree-level upper bound on the Higgs mass. 



m, 



< M| (cos^ 2(3 + 2.8 sin^ 2/3) , (6.9) 



generating more flexibility with respect to the experimental data. For example, for tan /? = 
2(4) one obtains ruh < 1.47(1.18) x Mz ~ 134(107) GeV. Using the above mentioned 



-20- 



possibility of 10 TeV for the high-energy, scale ||5^, producing < (1.96)^, the result 
would be ml < M| (cos^ 2/? + 13.2 sin^2/3). In this case, for tan/3 = 2(4) one obtains 
nih < 2.96(1.92) x Mz ~ 270(175) GeV. 

Concerning the rest of the spectrum, the charged Higgses are mixed with the charged 
sleptons generating the 8x8 charged scalar mass matrix written in eq. (A. 27). Nevertheless, 
similarly to the neutral scalar mass matrices where some sectors are decoupled, the 2x2 
charged Higgs submatrix is decoupled from the 6x6 charged slepton submatrix. 

The neutralinos are mixed with the right- and left-handed neutrinos producing the 
10x10 neutral fermion mass matrix written in eq. ( |A.49 ). As discussed in Section |5|, 
three eigenvalues are very small corresponding to the neutrino masses. The other seven 
eigenvalues arise from the mixing of neutralinos and right-handed neutrinos. 

As discussed also in Section ^, although the charginos mix with the charged leptons 
giving rise to the 5x5 charged fermion mass matrix written in eq. ( |A.47 ), the 2x2 chargino 
submatrix is basically decoupled from the 3x3 charged lepton submatrix. The former is 
like the one of the MSSM provided that one uses fi = Xivf- 

Let us finally mention that the squark mass matrices are written in eq. ( A.41| ). When 
compared to the MSSM case, their structure is essentially unaffected, provided that one 
uses /U = Aji^f, and neglects the terms proportional to Yy. 

For a more detail discussion of the characteristics of the spectrum we need more in- 
formation about the parameter space. As an example, let us consider the viable region 
studied in Fig. ^ with A = 0.1 and k = 0.4. We show first in Fig. ^ the masses of 
the CP-even neutral scalars as a function of the right-handed sneutrino VEVs. For this 
parameter space we can see from Appendix |A.1.1| that the mixing between the Higgses and 
the right-handed sneutrinos is of the order of a\^Vu = AxXvu, and therefore small compared 
with the relevant diagonal terms AjAjU^f^ = 9A^?;'^'^. Thus we have essentially doublet-like 
Higgses and the LEP bound for the lightest Higgs mass applies. The masses of the heavy 
and light Higgses, H and h, are shown in the figure with green dashed and solid lines, re- 
spectively. Concerning the former, its mass varies between 1748 and 2935 GeV. Concerning 
the latter, since tan/3 = 20 the upper bound is like in the MSSM, as discussed above. For 
the values of the parameters used in this example, we obtain m/j ~ 115.5 GeV. If instead of 
At = 1 TeV, we would have consider the 'maximal mixing' scenario |53], which in our case 
is obtained for At ~ 2.4 TeV, we would have obtained rrih ~ 126 GeV. As discussed also in 
eq. ( |6.8D , larger values can be obtained avoiding as much as possible the small mixing of 
the light Higgs h with the right-handed sneutrinos. For example, for A = 0.05 one obtains 
rrih ~ 117.5 GeV. Imposing in addition the maximal mixing scenario, rrih ~ 128 GeV. 

The three right-handed sneutrinos are essentially degenerated (up to small contribu- 
tions due to neutrino Yukawas), and we show their masses with a black dashed line which 
varies approximately between 357 and 1346 GeV. Let us remark that in general to obtain 
singlet-like Higgses, thus scaping detection and being in agreement with accelerator data, is 
also possible for small values of k. This can be qualitatively understood from the expresion 
of the corresponding mass matrix. In particular, the terms ^"^-pc-^R^yc-^R are of the order of 
K?v'^'^, and become very small when k decreases. 
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Figure 8: Masses of the CP-even neutral scalars as a function of the right-handed sneutrino VEVs, for 
the parameter space of Fig. 7b with A = 0.1 and k = 0.4. The gray and violet areas are excluded by 
the existence of false minima and tachyons, respectively. 
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Figure 9: The same as in Fig. 8 but for the masses of the neutral fermions. 
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Concerning the left-handed sneutrinos Vi in Fig. we see in the Appendix that their 
masses are basically determined by the corresponding soft masses, mj . Notice that the 
other terms in M^~.j^^~ -^^ are proportional to Vi or Yy^., and therefore negligible. On 
the other hand, the values of m~ are fixed by the minimisation conditions ( p. 61 ), and as 
a consequence they are essentially proportional to {Yu-/vi)v'^ for the viable region of the 
parameter space studied here. For example, for = 1 TeV in the figure, the values of the 
Yukawa couphngs are given by Yi,^ = 1.64 x 10"'^, Yy^ = 5.43 x 10"'^ and Y^.^^ = 9.85 x 10"'^. 
Using the VEVs Vi discussed above eq. (|6.1|), one obtains from the previous formula 

~ 1.8m,>j, and m^.^ ~ O.TTmp^. This can be checked with the figure. 

Let us finally remark that for the region of the parameter space discussed here, to work 
with other values of tan fj would not modify the spectrum obtained, with the exception 
of the masses of h and H. This is also true for the rest of the spectrum discussed below. 
For example, for tan /? = 5 we obtain essentially the same spectrum but with mn varying 
approximately between 1310 and 2332 GeV, and m/i ~ 112 (124 GeV for maximal mixing). 



It is straightforward to see from Appendix A. 1.2 that the masses of the CP-odd neutral 
scalars are very similar to those of the CP-even neutral scalars discussed above. In partic- 
ular, the masses of the pseudoscalar and left-handed sneutrinos are similar to the masses of 
the heavy Higgs and left-handed sneutrinos in Fig. The only differences appear for 
the right-handed sneutrino masses. Note e.g. that the terms 2ai^..^v'j^ and 2KijkKimk''^^ 



have different signs in eqs. ( A. 11 ) and ( A. 24 ), implying that now the masses vary approx- 
imately between 1 and 1.5 TeV. 



Concerning the charged scalars, we can see in Appendix A. 1.5 that the mass of the 
charged Higgs is very similar to the ones of the pseudoscalar and heavy Higgs. As mentioned 
above, the right- and left-handed charged sleptons are decoupled from the charged Higgs. In 
the Appendix we see that their masses are essentially determined by the corresponding soft 
masses, mg-c, mr . Although the former are free at the electroweak scale in our computation, 
the latter are fixed by the minimization conditions ( p. 61 ), and therefore we obtain the same 
masses as for the left-handed sneutrinos. 

In Fig. P we show the seven eigenvalues corresponding to the mixing of neutralinos and 
right-handed neutrinos. As mentioned in the previous Subsection, we have taken values 
for the soft gaugino masses that mimic at low scale the results from a hypothetical unified 
value at the GUT scale. In particular, we have assumed M2 = 1 TeV and consequently 
Ml ~ 500 GeV. As we can see in the figure, and can be deduced from the matrix ( |A.51 ), 
for the values of the parameters analysed we obtain almost pure Wino, Bino, Higgsino 
and right-handed neutrino states. The blue dashed (solid) line corresponds to the Wino 
(Bino) mass, which is determined approximately by the soft mass M2 (Mi). The Higgsino 
masses are determined approximately by the effective // term, Aji/? = 3Az^'^. We show with 
a green dashed (dot-dashed) line the heaviest (lightest) Higssino H2 {Hi). Their masses 
vary between 267 (242) and 617 (464) GeV. Finally, the three right-handed neutrinos 
are degenerated with a mass that can be approximated as 2ki''^. This is shown with a 
black dot-dashed line in the figure varying between 686 and 1620 GeV. Although in the 
present case the lightest neutralino is a Higgsino, due to our choice of input values with 
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Ml > SXv^, this can easily be modified by choosing other values of the parameters. The 
lightest neutralino can also be essentially a right-handed neutrino for small k. Let us 
finally remark that varying the values of the parameters also the mixing of states can be 
augmented. This can be obtained by making the diagonal entries similar to each other 
and/or increasing the off diagonal entries. 

On the other hand, from the 2x2 chargino submatrix in eq. ( |A.48 ) we can easily 
deduced that the mass of the charged Wino is approximately given by M2, and the mass 
of the charged Higgsino by the effective // term, jj, = Xivf. 

Finally, the eigenvalues of the squark mass matrices depend on the soft masses. As for 
the right-handed sleptons, in our computation these are free parameters at the electroweak 
scale. 



7. Conclusions and outlook 



We have performed the first detailed analysis of the /iz^SSM. As explained in the Intro- 
duction, this model was proposed |^5| as a SUSY standard model for solving the crucial 
/X problem of SUSY constructions, generating at the same time the small neutrino masses 
through a dynamical see-saw at the electroweak scale. This is due to the inclusion of three 
generations of right-handed neutrino superfields and the corresponding new gauge invari- 
ant couplings, eab^i uf H°[H^ and Kijki'fi'ji'^. The latter couplings break R-parity explicitly 
and therefore the phenomenology of the //z^SSM is very peculiar and different from other 
models, not only from those conserving R-parity, but also from those were R-parity is also 
broken. 

In this work we have extended the analysis of ref. |^5[, where the characteristics 
of the /xz^SSM were only introduced, and several approximations were considered in the 
phenomenological discussion. In particular, only one generation of sneutrinos were assumed 
to acquire VEVs. Here we have worked with the full three generations. We have written 
for the first time the corresponding scalar potential and minimized it in order to study 
the electroweak symmetry breaking. One-loop corrections have been taken into account in 
the computation. In total eight fields acquire VEVs. They are, in addition to the usual 
Higgses, the right- and left-handed sneutrinos. Notice that minima with some or all of the 
VEVs vanishing are in principle possible, and therefore one has to check that the minimum 
with non-vanishing VEVs breaking the electroweak symmetry, and generating the fi term 
and neutrino masses spontaneously, is the global one. 

Obviously, due to the many VEVs and the new couplings, the parameter space of 
fivSSM is very involved. After discussing in detail the strategy to follow in the low-energy 
analysis, we have studied viable regions of the parameter space which are left after imposing 
several constraints. In addition to discard regions with the false minima mentioned above, 
we have discarded also regions with tachyons, as well as those where the Landau pole 
constraint on the couplings at the GUT scale is not fulfiled. Of course, reproducing neutrino 
data is also used as a constraint in the parameter space. Results are shown in Figs. 2-7. 

Finally, we have discussed the particle spectrum. The breaking of -R-parity generates 
complicated mass matrices and mass eigenstates. The presence of right and left-handed 
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sneutrino VEVs leads to mixing of the neutralinos with the neutrinos producing a 10x10 
matrix. Indeed three eigenvalues of this matrix are very small, reproducing the experimen- 
tal results on neutrino masses. On the other hand, the charginos mix with the charged 
leptons giving rise to a 5x5 matrix. Nevertheless, there will always be three light eigenval- 
ues corresponding to the electron, muon and tau. Concerning the scalar mass matrices, the 
neutral Higgses are mixed with the sneutrinos, and the charged Higgses with the charged 
sleptons, and we are left with fifteen (eight CP-even and seven CP-odd) neutral states and 
seven charged states. Notice however that the three left handed sneutrinos are basically 
decoupled from the Higgs-right handed sneutrinos, and also the six charged sleptons are 
decoupled from the charged Higgses. 

Given the interest of the lightest Higgs boson mass in the analysis of SUSY models, 
we have discussed in detail the mass of the lightest CP-even neutral scalar in our model. 
The upper bound turns out to be similar to the one of the NMSSM, about 140 GeV after 
imposing the Landau pole constraint up to the GUT scale. For the precise masses of the 
Higges and of the rest of the spectrum, it is not possible to give a result valid for the 
whole parameter space, given the complicated structure of the model. Nevertheless, we 
have pointed out several interesting characteristics, and analysed particular regions and 
possible variations. An example of a possible spectrum is shown in Figs. 8-9. 

Once we have checked explicitly that the parameter space of our model contains viable 
solutions and the associated spectrum is interesting, and given the hope that the LHC will 
be able to test SUSY, it is then important to study in detail the collider phenomenology 
of the /Ltz/SSM. In particular, the impact of the new couplings on the usual SUSY searches, 
and indeed novel signals that might facilitate the confirmation of the /xi/SSM as the ad- 
equate SUSY Standard Model. This necessary task will be the subject of a forthcoming 
publication. 
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A. Mass matrices 

In this Appendix we will study the general mass matrices generated in the /UZ^SSM. For 
this study we will use the indices k,l,m = 1, 2, 3, and a, /3, 7, 5 = 1, 8. 
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A.l Scalar mass matrices 



Here we study the scalar mass matrices. Let us recall that concerning the Higgses, the 
neutral ones are mixed with the sneutrinos, and the charged ones with the charged sleptons. 

A. 1.1 CP-even neutral scalars 

The quadratic potential includes 

^quadratic = ^'a-M^^^S'/j + ... , (A.l) 

where = {ha, hu, (^f)^, (i^i)^) is in the unrotated basis, and below we give the expres- 
sions for the independent coefficients of Mg^^ 

Mlha = "'H, + -^{^Vd -^1 + ^rn} + Wjyly] + Wvl , (A.2) 

-2Y^..XjVdiyi + Y^.^Y^..u^i^^ + Y^.^Y^j^iyiiyj , (A.3) 

^hd{v9)R = -aXiVu + 2XiXjVdiyj - 2XkKijkVuiyj - Y^..XkVjU^ - Y^.^XiUjU^ , (A.5) 

= -(^Xi'"d + auji^j + 2XiXjVui^j - 2XkKiikVdi'i + 2Y^.^KiikVjuf + 2Yy.^Y^^^VuVk > 

(A.6) 

= \G^^d^r - {Y.^^^jvl + Y.^^Xkulv'^) , (A.7) 
^huiyi)^ = ^''ij^j ~ -Y^nn - 2Y^^.XjVdVu + Y^.^KijkViy'^j + 2Yy^.Y^^^VuUk , (A.8) 

q2 2 

M|.)H(i;.)« = + —UiUj + -G\uki^k + vl- vl)dij + Y^.^Y^.^vl + y,.,F,.,i/^z/f , 

(A.9) 

M(2~ = a„,.Vu - Y^.-Xkydi^k - ynk>^jVdT^k + "^Ynui^jikVuT^i 
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^lv'r)R{V'^)R = n^lf. + ^O'Kijk^k ~ '^^kl^ijkVdVu + '^l^ijkl^lmk^lt^m + '^l^ilkl^jruk^l^m 



+\iXj{vl + vl) + 2Y^^^KijkVuJ^i - {Y^^^^Xi + Y^^.Xj)vdVk + Y^^.Y^^.vl + Y^^.Y^^.i^ki^i . 

(A.ll) 

Then the mass eingenvectors are 

Sa = KpS'p , (A. 12) 

with the diagonal mass matrix 

{Mflff = K^Ml^R^s ■ (A.13) 
A. 1.2 CP-odd neutral scalars 

In the unrotated basis P'^ = {Pd-,Puj (yiY-, (^iY) we have 

Vquadratic = P'aM|.^^ P'^ + • ■ • (A. 14) 

Below we give the expressions for the independent cofHcients of Mj, ^ 

-Tf,,.XpS^i + Y.^^Y.^.vlv'l + Y^^^Y^.^ViV^ , (A.16) 
M|,^P^ = ax,i^f + XkKijki^fu^ , (A. 17) 

Mp^(pcy = aXiVu - 2XkKijkVui^j - Y^,.^XkVk^j + Y^.^XiU^Vj , (A.18) 

Mp^(~e)7 =aXiVd - a,,.Mj - IXkKiikVdi'i + 2Y^,.^KiikVjyi , (A.20) 
^PuiJ^iY = -""'ij^j - ^ikf^ljk^l^'j ' (A-21) 

^h'{-,y = ""i, + \G\^kyk + vl- vl)5ij + Y,^^Y,^y^ + y.,,y.,,i/^i^f , (A.22) 

(A.23) 
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+ h>^A'"d + vl) - (i^i^fciAj + Y^^.\i)vdyk - 2Y^^^KijkVuJ^i + Y^^.Y^^.vl + Y^^^.Y^^MkUi . 

(A.24) 

Then the mass eingenvectors are 

P„ = <^PV , (A.25) 

with the diagonal mass matrix 

(^pj)' = <Ks^PS ■ (A.26) 

A. 1.3 Charged scalars 

We give here the mass matrix coefficients for the charged scalars which follows from the 
quadratic term in the potential 

^quadratic = S'^M^i 8;,+ . (A.27) 

a.j3 

The unrotated charged scalars are S'^ = (i?j',i7+,e^,/I^,T^,eJ,/U^,T^), and 
1 

MH^Ha = ^la + 2^2^^ - I'ii'i) + -ji^i^i +^1- vl) + AiAjZ/>; + Y^.J^.^ViVj (A.28) 
1 

^l^H^ = ra\^ + 2^2(^1 + ^i^i) - + - <) + A^A,I/fl/^^ + Y,^Y,^^v]vl (A.29) 



^^d^^.. = '^^i^i' + 2^2VdVu - k><iVdVu + hKijkffi^j + Yi,..\jVuiyi (A.30) 

+ y.,y.^,z.,v^ + ye«i;,,^i (A.31) 

= «e».^<^ - 5:^e,Afci;„z.^ (A.32) 

^«,e., = Ml,e,^ (A.33) 

^i,eH. = + fl-l^fei^fc - + + ye,,i;,,^^d + ^H^e^.^feZ^i (A.34) 

Ml^H, = ^^VdH - Y,,.XkU-,u^ - Y,,.Y,^.VdUk (A.35) 
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Ml^^H, = -ae,,yj - Ye,Y,^.VuU^ (A.37) 

Ml^H^ = -YeJXjUkU^ + Y^^.vau^) , (A.38) 
where Oe^^ = {A^Ygjij. Then the mass eigenvectors are 

= <S'J , (A.39) 

with the diagonal mass matrix 

{Mfr'fap = Rt,Ml^ Rfs ■ (A.40) 

It is worth noticing here that if we allow the presence of the lepton number violat- 
ing terms in the superpotential, Xij^LiLje'^, discussed in the Introduction, they would 
contribute to the above charged scalar masses. 

A. 1.4 Squarks 

In the unrotated basis, u'i = (^^^,2^.) and d'i = {dLi,d*j^), we get 

Kiuadratic = u' + \d'^ Mj d' , (A.41) 

where 

/ M| M| 

M| = '^^^^ I , (A.42) 

with q= {u',d'). The blocks are different for up and down quarks, and we have 



and 



= (^ui^'^u - YuijXkVdi^k + YuikYui^i^ii^k , 



dL^Lj Qij 6' 2 2 
M| =w?~ , (A.44) 
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where a„^^- = {AuYu)ij and a^^. 



with the diagonal mass matrix 



{AdYd)^j. 



For the mass state we have 



diag\2 



(A.45) 



(A.46) 



It is worth noticing here that if we allow the presence of the baryon number vio- 
lating terms in the superpotential discussed in the Introduction, ^ijk^iQjd%, they would 
contribute to the above squark masses. Actually, even if they are set to zero, one-loop 
corrections will generate them, as discussed in Appendix However, these contributions 
are negligible. 



A. 2 Charged fermion mass matrix 



Charginos mix with the charged leptons and therefore in a basis where ^ 

" - 



and ^ = {—iX ,ej^, ^j^,Tr ), one obtains the matrix 



where 



Mc 



92V u 



A/J 
Mc 












\ 



Ye^^Vd Ye^^Vd Ye 



ei^Vd 
Ye2iVd 

Ye^-^Vd Y^z^Vd Ye,^^Vd j 



Ye2\Vd Ye„„Vd 



(-iA+,^+,e^,/x+,T+) 



(A.47) 



(A.48) 



/ M2 

92Vd 

92^2 

A. 3 Neutral fermion mass matrix 

Neutralinos mix with the neutrinos and therefore in a basis where = {B'^-, Hd-, Hu, fRi-, ^Li), 
one obtains the following neutral fermion mass terms in the Lagrangian 

1 



(A.49) 



where 



with 



M = 



Mr 











-Avd 


Avu 











\ 







A/2 




-Bvu 

















Bvd 







-Aiu„ 


-A2W« 








Avu 









-XiVd + Y^.^Vi 


-A2Ud + Y^y^^l^, 


-X-iVd + Y^^^v, 














-XiVd + Yy^^Vi 




















-X2Vd + Yy^^Vi 






2K23jVj 




V 








-AsWu 


-AsWd + Y^.^Vi 




2K32jl^j 




/ 



(A.50) 



(A.51) 
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where A = ^ sin 9w, ^ ~ ^w, and 



m 



/ ~Vl^^ vl^^ ^ "^''1*^''' ^Vll^^W ^V12^^U Yui3Vu\ 

V -Jl'^s :J|'^3 y^g,!^? y^3i?;„ y^32t;„ Y'i.ag?;^/ 



(A.52) 



B. Couplings 



In this Appendix we show the relevant couphngs involved in the computation of the one- 
loop radiative corrections to the scalar potential tadpoles and the CP-even scalars masses. 

B.l Scalar up squarks— up squarks 
With the definition 



r - „S'°u'u" oiO ~i ~i* J 



(B.l) 



we get 



.,3'°u'u" „S">u'u'- 
\iiaRiLj ilaRiRj 



(B.2) 



where 

S'°u'u'* 
iJaLiLj 

S'Ou'u'* 
9aLiRj 

„S'°u'u"- 
iJaRiLj 



up Sap (- 5 5^ + i 5'^) - 2 6i2 Vu Yu-i Yu^^ , 

- {AuYu)ij + 6al A; Y^^ - da-2,l Yui^l^m Yu,j + Sa-5,l {vdk " Yu^.l^m) Y^^ 



9aL~R., 



yaRiRj 



and we have defined 

up = {vd, Vu, V2, V3, vi, v^, 1^3) ; ^ij = diag(+, -, +, +, +, 0, 0, 0) 
while 5ij is equal to one for i = j, and zero for i ^ j. 

B.2 Scalar— down squarks— down squarks 

With the definition 



(B.3) 



(B.4) 



r _ ^S'Od'd'* <j/0 J' J/* , . . . 

iiaij '-'a "i ' ' 



(B.5) 



we get 



S'Od'd" 

Haij 



„S'°d'd'* „S'°d'd'*' 

iJaLiLj iJaLiRj 

„S'°d'd'* „S'°d'd'* 

\yaRiLj YaRiRj 



(B.6) 
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where 



= -'^"1 i^dYd^j + 1^1 Xi Yd,j + Sa-5,i vu Yd,. , 



S'°d'd'* ^ S'°d'd'* 



9aRX = -3 5" - 2 Sal ^d^d^ld,, • (B.7) 

We find the couphngs in the squark 5i,2 basis via fffij^^ = ^uidahl^ )-^Jm- 



B.3 Scalar— quark— quark 

With the definition 



^ = S'a Ui Uj + S'^ d,dj + -- - , (B. 



we get 



and 



Jaij "az Uij 



-6a2Yu,, , (B.9) 



9!^ = -SaiYd,, . (B.IO) 



B.4 Scalar— scalar— up scalars-up scalars 
With the definition 



^ - 9a3ij ^a^B^iUj H > (B.ilj 



we get 



\9al3RiLj 9al3RiRj 



where 



„S'°S'°u'u'* 
9al3LiLj 


= 5.p{-\9' + h9"] 


) - Sa2 Sp2 


Y Y 


9 a/3 Li Rj 




- Sa-2,l 5/3 


-5,m Yl,^^ Yu,. 


S'OS'°u'Z'* 
9apRiLj 


— 9al3LjRi ) 






9al3RiRj 


= -5^0/35'^ - Sa2 Sf32 


Y Y 





(B.13) 



-32- 




B.5 Scalar scalar down scalars down scalars 

With the definition 

^ - 9a0ii ^a^B^id^ H , (i^-14j 



we get 



/ S">S">d'd'* S'°S">d'd'* 
1 s'Os'Od'd'* „S'°S'°d'd'* ' 



where 



^S'°S'°d'd'* _ ? /' 1 „2 , 1 /2^i A A V V 



R'°S'°d'd- -A A \ V 

9al3LiRj - 2°°2 Of3-5,lM Yd.j , 

S'^S'^d'S" _ S'°d'd'* 

9al3RiLj — 9al3LjRi ) 

sS?"'* = l^^^P 9" - 6ai h2 Yd^^Yd,, . (B.16) 



1 



C. Tadpoles 



In this Appendix we present the leading one-loop DR tadpoles (i.e. the ones involving 
s(quarks) in the loop) which enter into the minimization of the neutral scalar potential 
(see Figirg), 



167r2 

where X = {u, d, u, d), and 

3 

^S'o =Y.3g'2^4mf^A,{ml) , (C.2) 



k=l 



4o = -Y.^9'jFAoiml), (C.3) 



k=l 



where f = u,d and is the 1-point Passarino-Veltman function [47|. 
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Figure 11: Self-energy diagrams 



D. One loop self-energies 

Here we list the leading one-loop DR self-energies of the CP-even scalar mass matrix 
represented in Fig. (|Tll), 

3 _ 2 

f=u,d k=l 

f=u,d k,l=l 
6 

f=u,d k,l=l 

where Ni is the number of colours, which is 3 for a (s)quark and Bq is the 2-point Passarino- 
Veltman function p^]. 

E. Re normalisation group equations of Yukawa couplings 

In this Appendix we give the RGEs of Yukawa couplings including Aj and Kjofc . Defining 



7^ = -'2{KukKjik + XiXj + Y^^Y^^.) , (E.l) 

7h: = lal + ^gl - 31;,, - XiX^ - Y.^.Yu,, , (E.2) 

ill = \9l + - - sy^^.y^,, - a, a, , (e.3) 

= + - Ye^^Ye,, - Y,^,Y,^, , (E.4) 

7fr=^H; = -^^.A, , (E.5) 

^% = \g\-2Y,^,Y,^, , (E.6) 
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Figure 12: One-loop generation of the X'-j^LfQ'jd'j, term in the superpotential. Note that it is propor- 
tional to Fy, Yd, and A. 



at one- loop level we have the following RGEs: 



dt 

^Y 
dt 



d 






(E.IO) 


1 

' ~ 167r 






(E.ll) 


1 

167r2 






(E.12) 


^Y 

dt 






(E.13) 


-Yd - 
dt 






(E.14) 


Ay - 






(E.15) 



where t = — InQ, with Q the renormalization scale. 

It is worth noticing here that one-loop contributions in the //i^SSM will generate one 
of the usual lepton number violating terms mentioned in the introduction, A^^-^L^Q^fi^, as 
shown in Fig. The corresponding RGEs are: 



However, this contribution is proportional to the neutrino Yukawa coupling, and therefore 
can be neglected in the computation. 
Finally, for the VEVs we have 

1 d 



1 d 



167r2 dt 



Y^|-. = -.7£-^^.7fr ' (E.19) 
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